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Abstract 



We extend the model of Ancilla Driven Quantum Computation (ADQC) by considering gates 
with arbitrary entangling power. By giving up stepwise determinism, universal QC can still be 
achieved through a variable length sequence of single qubit gates and probabilistic "repeat-until- 
success" entangling operations. This opens up a new range of possible physical implementations 
as well as shedding light on the sets of resources sufficient for universal QC. 
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4^ , 1 Introduction 



A fundamental question is "what set of resources are required for universal quantum com- 
putation?". Many models have been proposed ranging from the conventional gate-based [1], 
^ measurement-based [2], adiabatic [3 , and topological [4 quantum computation. These mod- 

els utilise different sets of resources and are suited to various physical implementations. By 
finding new models with different sets of sufficient resources, this may open up new ways 
of implementing quantum computation as well as enrich our understanding of quantum 
computation itself. 

\^ A hybrid model combining aspects of the gate-based and measurement-based approaches 

was introduced as Ancilla Driven Quantum Computation [5J. The key feature of ADQC is 
the restriction of resources to a single unitary interaction, and direct access (initialisation and 

CO measurement) only to an ancilla qubit that can be coupled sequentially to system (register) 

qubits. ADQC uses entanglement between the ancilla and register, and the kickback induced 
by measurement on the ancilla to drive unitary evolution of the register. By coupling 
^ the ancilla to various register qubits and choosing different measurements on the ancilla, 

universal QC can be achieved. 

rN There is often a trade off between easy access and manipulation or long coherence times. 

^ ADQC lends itself to physical systems where register qubits with long decoherence times 

are difficult to manipulate while relatively short-lived ancillary systems are more easily 
controlled and can be quickly initialised and measured. There has been much work that looks 
at dealing with such properties with the use of ancillas for specific physical implementations. 
Work on optical clocks using aluminium ions have employed magnesium or beryllium ions 
as an ancillary system to account for a lack of optical accessibility with aluminium ions |6]. 
Ohshima [7 considered maintaining low decoherence of quantum dots by only activating 
access through an ancilla qubit in the same cell. It has been proposed that isolated, stable 
NV centre nuclear spins be used as qubits manipulated by neighbouring electron spins [H [9] 
and Bermudez et al developed an proposal for nuclear spin interactions mediated by electron 
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spins that effects an Ising type interaction [9 . Ion trap+photon systems such as in [10 and 
sohd state systems with bahistic electrons have been considered for a class of systems that 
involve generating quantum gates through scattering between a static and flying qubit [TT]. 
In such cases we may have no interaction-time tuning [11 or restricted range of dynamic 
modulation which is particularly relevant for the work of this paper. 

These proposals often look at one system or particular parameters with a focus on gate 
based QC. The proposals of Anders et al [5 and Kashefl et al [12 from which we will develop 
our proposal considers a scheme without reference to a speciflc physical system. A general 
description is provided, of the minimal resources require to achieve universality, expressed 
as a flnite resource set that can be compared to Measurement Based Quantum Computation 
(MBQC) and the resources to form cluster states but aimed towards a Gate Based Quantum 
Computation (GBQC) style circuit. 

In the original ADQC scheme, the computation proceeds step- wise deterministically up to 
Pauli corrections on the flnal output, very much in the spirit of MBQC [2[[13]. This requires 
that the entangling unitary between register and ancilla be of either of two speciflc gates, 
H ^ H.CZ or SWAP.cz. No other entangling gate would allow the measurement-induced 
kickback to be unitary, and permit the different computation branches corresponding to the 
possible measurement results to be reunited via Pauli corrections [5 . 

Here, we show that by relaxing the requirement for stepwise determinism a much larger 
set of entangling gates can achieve universal QC which may open up a range of physical 
implementations with flxed arbitrary coupling strengths and interaction times. Investigation 
of to what extent the required entangling properties can be relaxed has been performed for 
MBQC [14, J^. Unlike that, we do not use many-body physics techniques but expand the 
resource set using GBQC concepts. 

Pauli corrections and the corresponding direct access of the register will not be required. 
This may also reduce the control requirements for characterisation as described in [16 . This 
is achieved by gate approximation, repeat-until-success strategies, and multi-step measure- 
ment. Single qubit unitary gates can still be implemented deterministically while two qubit 
gates and measurement and state initialisation will be probabilistic and require the devel- 
opment of probabilistic protocols of which we provide examples. 

^^2 Overview of ADQC 

We will review how ADQC implements single qubit unitary gates and which conditions 
on the available resources are necessary for the method to be step- wise deterministic. We 
highlight the resource requirements that our proposal will extend; a full description can be 
found in EHHHT]. 

The evolution of a qubit in a quantum register can be driven by preparing an ancilla 
qubit system in a state |a), coupling that ancilla with an entangling interaction. Ear and 
then measuring the ancilla in some basis {|m+), |m_)}. After measurement, the evolution 
of the register qubit can be described by a Kraus operator K± = {m±\E\a) \lS\ [T9l [20]. 

The interaction Ear can generally be composed of a product of gates that act locally on 
the individual systems and a gate that produces some entanglement between the systems. 
For classiflcation of two qubit gates, we turn to the canonical decomposition [211 [221 123] : ^ 
unitary that acts on a system of a pair of qubits A and R, Ear^ can be expressed as 



Ear = {Va VR)A{a 



(1) 
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where Ua^ Ur^ Va, Vr are unitary gates that act only on the local systems of qubits 1 and 2 
and A{ax, ay, a^) = exp(-z(aa,cra, cr^; + ayCFy + a^a^ (g) a^)). 

If the final action on the register is to be unitary then the Kraus operator must be 
proportional to unitary such that K±K\. = p±l where p± is the probability of the meas- 
urement. This does not depend on the local unitary gates acting on the register. The 
{ax^ocy^az) parameters are the only parameters unique to Ear that determine whether 
the measurement back-action will be unitary. The preparation of the state of the ancilla 
system can be adapted for any unitary that immediate follows it; Ua is effectively under 
full control and can be represented entirely by ancilla state preparation. Similarly choice of 
measurement basis equates to freedom of choice of the unitary after interaction, Va- Only 
the /S.{ax^ay^az) component is of interest when considering the entangling capabilities of 
the interaction [2X1122]. 

K± are generated probabilistically and are not equivalent. In order to be deterministic, 
a key idea from MBQC is utilised: if PK- = Kj^ where P is a Pauli operator correction, 
this correction can then be commuted through several applications of the Kraus operator 
and local unitary gates. Different computation branches associated with the measurement 
outcomes can be reunited by a Pauli correction. 



{VRPK_UR).(yRPK_UR).(yRPK_UR) = P'P"P"\VrK_Ur).{VrK_Ur).{VrK_Ur) 



To have the capability to enact any arbitrary unitary gate on the register, a sequence of 
VrK-Ur should be universal for single qubit unitary gates. 

2.1 Control-Z Hadamard example 



H Figure 1 Anders et al. depiction of ancilla-driven implementation of a single qubit rotation [5]. 
The ancilla and register qubits are coupled with CZ and the local unitary gates are chosen such that 
the interaction remains symmetrical with respect to ancilla-register exchange. A rotation J{/3) 
is enacted on the register a result J{/3) is enacted on the ancilla which is then measured in the z 
basis with a result j=0,l. 

An example can be performed with Ear = {Ha Hr).CZar- The ancilla is prepared 
in the |+) state, couples with the register qubit with Ear and then undergoes a unitary 
J(/3) = HR^{j3) before being measured in the computational basis. The action on the 
register qubit for a measurement result \j) is Xj^Jr{I3). The difference between the two 
resulting unitary operations can be corrected by Xr. Any arbitrary single qubit unitary can 
be decomposed into four rotations e*" J(0) J(/3) 7(7) J(^) so to implement any arbitrary single 
qubit unitary up to a global phase, four ancilla interact ion- measurements are performed 
changing the parameter of the local unitary on the ancilla to be the Euler angles (5, 7, /3 and 
then 0, 




CZ,,i, 



X\J{0)X^J{l3)X''J{j)X^J{S). 



(2) 



4 



Ancilla Driven Quantum Computation with arbitrary entangling strength 



As in MBQC, the Pauli corrections can be commuted through each apphcation of J(/3) if we 
make adaptations to the local unitary applied to the ancilla measurement basis [21 [131 [12 • 

X'j{0)X^J{/3)X^J{-f)X^J{S) = X'Z^X^Z^J{0)J{{-l)^/3)J{{-iyj)J{S) (3) 

To extend to larger computations all corrections on register qubits are required to inter- 
change with future entangling operations A{ax^ ay^az) so that they remain local corrections 
on the register and on the ancilla preparation or measurement basis choice. This allows all 
the Pauli corrections to accumulate at the final step. This provides the condition that the 
entangling operation A{ax^ay^ az) tensor commutes with the corrections [5 ; as a result only 
two classes of coupling are universal. 

2.2 Conditions on the class of interactions 

In order to allow achronical Pauli operator corrections ensuring general stepwise determ- 
inistic single qubit gates exactly, we must have interactions that belong to the classes of 
interactions that are locally equivalent to CZ or CZ -\- SWAP gates. Expressed in terms 
of A{ax^ay,az)^ these classes are A(|,0,0) and A(J, |,0). There exists broader classes of 
interactions which fulfil the weaker condition that the Kraus operators acting on the register 
qubit are proportional to unitary. These are the classes for which there is at least one 
parameter equal to zero- A(a, 0, 0) and A(ai, 0^2, 0) up to symmetries. It is symmetrical up 
to local unitary gate corrections with respect to ± exchanges, | shifts and refiections in J 
acting on the parameters {a ) [21] ^ . The values are also symmetrical with respect 

to permutations [23 . This allows us to consider only cases where \aj \ < J and to classify an 
interaction by the number of non-zero parameters, e.g. e-*("crz(g)cr^) equivalent, up to local 
unitary corrections to any case where (a^ = 0, aj = 0, a/e 7^ 0). We will demonstrate a way in 
which we can eschew a stepwise construction of unitary gates, thus not requiring Pauli cor- 
rections and allowing the broader class A(a,0,0), locally equivalent to the Control-unitary 
set of gates, for ancilla driven quantum computation with arbitrary interactions. As a cost, 
we will not be performing exact unitary gates but efficiently generated approximation. 

Some characteristics are general for any member of the class so we will use e~'^^^^^^^ 
to represent them. Other effects are dependent on specific interactions with descriptions of 
the local unitary gates. In the next section we will explain our choice of interaction in those 
cases. 

3 Single qubit gates using {H (g) H).C-T 

An interaction A(a,0,0) where < a < J has a restricted set of conditions on the choice 
of ancilla preparation state \a) and measurement basis that allow for a measurement 

induced back action represented by the Kraus operator (m^ | A|a) to be unitary. The resulting 
set of unitary gates that can be implemented deterministically is smaller than for A( J, 0, 0). 
Using A(a, 0, 0), a < J, only two gates can be implemented deterministically. If we consider 
just using the non-local part of the Cartan decomposition e"*"^^*^^^, the only two possible 
unitary gates that can be enacted independently of the random measurement result are 
achieve by preparing the ancilla in the computational basis. A |0), |1) ancilla input always 
corresponds to a Rz{2a)^Rz{—2a) unitary gate on the register qubit respectively. 

However, with local unitary gate corrections any member of the class A(a,0,0) is equi- 
valent to a Control-Unitary gate. Consider this in conjunction with a fixed local unitary 
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post-action, Ub^ on the register. Uh^R-C-Ua can implement a two gate finite set {Uq^Ui} 
where = t/^, Ui = UaUb. If the Lie algebra closure of Uq and Ui covers SU(2) then these 
two gates form a universal set for single qubit unitary gates 



Ua 



Ub - - U, 



X 



u - 



H Figure 2 A generalised description of the action by which members of the interaction class 
A(0 < a < f ,0,0) can generate two unitary gates whose commutation rules do not form a closed 
set. Uo = Ub and Ui = UbUa 

To demonstrate our proposal we choose the fundamental unitary for ancilla driven quantum 
computation E = [H H).C-Ri (^) (which corresponds to a — This specific gate 

is chosen for two reasons; 1) It is directly comparable with the gate Ear = {Ji ^ H).CZ 
from |5| but with a smaller rotation angle parameter of the Control-unitary, 2) It will 
generate {T, HT} which generates the same group as the finite set {H^ T} which is well 
studied and proven to be universal [23] and of interest for its applications in fault tolerant 
quantum computation [2H |25] . The method can then be generalised to arbitrary coupling 
strengths and local unitary gate products. With this method, the circuit is programmed by 
the ancilla state preparation and requires no further manipulation nor measurement of the 
ancilla. This may have benefits for particular physical systems depending on the lifetime 
and robustness of the ancilla. 

With this choice of interaction for producing single qubit gates, it is necessary to demon- 
strate the ability to perform measurements, initialise the register qubit into a specific state 
and enact a two qubit entangling gate in order to achieve universal quantum computation. 
We will address these issues in later sections. 



3.1 Two parameter interactions and stochastic ADQC 

The above method can only be employed with "one-parameter" interactions — the class 
A (a, 0, 0). There does not exist an ancilla preparation basis for which the resulting action is 
deterministic in the case of the A{ax^ o^y^ 0)- "two-parameter" class. Instead there must be a 
measurement of the ancilla which results in an action that is dependent on the measurement 
result. 

A two-parameter class interaction can be seen as a succession of one-parameter class 
interactions with local unitary corrections to account for permutation of the parameters 



e.g. 



- i ( Q! a: (Ja; (g) (Ja; + a z Cr^ ) 



H)C-Rz{4:ax) (up to local unitary corrections). Viewing it this way demonstrates why it 
should not have a strictly deterministic set of parameters but also why the form of the 
unitary actions brought about is always R2{/3)Rx{"y) = HJ{/3)J{^)H. Considering this the 
issue for using the two parameter class is just the inability to make a deterministic choice 
between two pairs of 7} to produce a universal single qubit gate set. It is similar in effect 
to flipping the resource requirements for the one-parameter case in time so that rather than 
preparing in the {\j)} basis, we only measure the ancilla in that basis. The resulting gates 
on the register would still be the same {Uq^Ui} as described before but randomly generated 
with probabilities po^Pi dependent on the initial ancilla state. 
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In fact, we can generate an approximation to any arbitrary single qubit unitary gate in 
this way. The stochasticahy generated sequence Ylk Ui{k) wih perform a random walk on the 
compact set of unitaries and be guaranteed by Poincare recurrence to reach an approximation 
of any unitary eventually. The generation of the strings Yl^ Ui(^k)^ though probabilistic, are 
still consistent with the conditions of the Solovay-Kitaev theorem [18] and so the efficiency 
of these approximations might be improved. 

Future research will look at finding the gate count required to hit a target averaged over 
all unitaries for a given error and the scaling with the error. In addition to the relevance 
to our work, it may serve as a benchmark for any circuit compilation sequence: to see how 
much better inputting it through a string of ancilla is compared with than "measuring them 
all and letting God sort it out". 

As a preliminary investigation, we simulated the random generation of gates using the 
interactions {H (g) i^)A(0,0, f^) and A(y^,0, y^), ancilla preparation state |+) and meas- 
urement in the computational state basis. For the one-parameter class interaction the 
resulting unitary gates are Uq = HR^{^),Ui = HR^{—^), po = pi = ^; for the two- 
parameter class, [/q = i?^(|)i?^(|),/7i = Po = Pi = \' The target unitary 
was Ut — Rx{^)' At each step a gate corresponding to the {Uq^Ui} of each interaction 
was multiplied to the product of the previous step starting with the identity operator. The 
normalised trace distance of the product, F = ]\k ^i{,k)i ^tnd the target unitary was evaluated 
at each step until within a chosen error size e < 0.05. 



\\V-Ut\ 



2-\Tx[VWt]\ 



< e 



(4) 



The number of gates required for this to occur was collected 1000 times and used to create 
a probability distribution for the gate count required to reach the target unitary (see figure 



probability density 




10000 20000 30000 40000 50 000 60 000 

(a) iJHA(0,0,f^) 



probability density 
IJ E-4 




10 000 

(b) A(^,0,fL) 



Figure 3 Probability distribution of required gate count to achieve target Ut = Rx{^)- a) Use 
of a single parameter interaction in a 20 bin histogram, b) Use of a two parameter interaction 
in a 20 bin histogram. The probability distribution corresponding to the exponential distribution 
parametrised by the mean of the results is displayed by the solid red curve (colour online). 



This small investigation reveals some behaviours to be considered in further research. 
While the target is not achievable in a single step and there is no finite probability for 
success per step, the aggregate behaviour over many steps, taken over a large number of 
simulations, can be modelled as a geometric or discretised exponential distribution. This is 
true for both one-parameter and two-parameter interactions. Particular target unitaries may 
cause anomalous effects; the two-parameter case is able to produce an exact solution of the 
target unitary in 4 steps which causes a large peak in the distribution and then suppresses 
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the probability of a result for several steps after (see figure 13d in Appendix [B]) but with a 
large enough bin size the exponential model dominates. This provides some features to test 
when extending the average over a large number of target unitaries. An average over many 
random target unitaries may smooth out such effects and verify the general applicability of 
the exponential distribution model. Under this model, we can calculate the probability of 
reaching any target unitary within a fixed number of steps and increase the number of steps 
until it provides a desired fidelity. 



L.4 Entangling gates 

Given the ability to implement any single qubit unitary gate, universal quantum compu- 
tation requires that we are also able to implement an entangling two-qubit unitary gate 
between register qubits. Note that direct interaction between register qubits not allowed, 
the interaction must be mediated by an ancilla qubit in the ADQC model. ADQC is capable 
of implementing an entangling gate by the use of one ancilla and two implementations of 
the same fundamental interaction Ear as used in the implementation of single qubit gates. 
The measurement of the ancilla results in two outcomes but they are equivalent up to local 
unitary corrections. 




-u- 



T - H - y' - \^)- 



















y 



H Figure 4 The circuit for implementing a two qubit gate and its generalisation. Vi V2 allows 
us to make any local unitary corrections to the register qubits. The preparation of the ancilla state 
is represented by Ui while the choice of measurement basis is represented by Ug. 

In our proposal, we attempt the same use of a single ancilla and two interactions but with 
weaker coupling strength A(0 < a < J). We have the freedom to prepare and measure the 
ancilla in any state/basis, as well as perform local unitary operations and post corrections 
on it and the register qubits. Due to this, we will focus on the e"*"^^'^^^ form of the 
interaction, with any local unitary corrections being incorporated into a single gate; the 
local H post-interaction unitary on the ancilla can be removed by the appropriate choice of 
Ua — see figure |4] . 

4.1 Interpreting the unitary and entangling conditions 

The operation on the register qubits after measuring the ancilla must be unitary and en- 
tangling. This will restrict the unitary operations, prepared states and measurement bases 
we can use on the ancilla. We need to be able to write the unitary and entangling conditions 
as some parameter restraints on the circuit. Because the interaction is of the class A (a, 0, 0), 
we can make some simplifications. If a register qubit is in a computational basis state 
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z = 0, 1, then it will, through the interaction A(a,0,0) cause a unitary — l)*2a) action 
on the ancilla. So for two register qubits there are four potential final ancilla states corres- 
ponding to the computational basis |aoo) for |00)i2 etc.. Therefore the transformation of a 
general register qubit and ancilla state, |a)|^)i2 will be transformed by the circuit thus: 

|a)|^)i2 = \a)J2Cij\ij) 12 ^ X]Qi|aii)Kj)i2 (5) 

Each final ancilla state is a unitary evolution of the original determined by the parameters 
of the circuit: 

\ai,) = UgR,{{-iy2a)UaRi{{-iy2a)\a) (6) 

After a measurement of the ancilla in a state |m) the (unnormalised) state of the register 
will be 

'^{m\aij)Cij\ij) 12 (7) 

The evolution of the register can be represented by a Kraus operator in a matrix represent- 
ation: 

Km = diag((m|aoo), (m|aoi), (m|aio), (m|aii)) (8) 

For this operator to be proportional to unitary must be the same Vz, j. This 

encapsulates an expression of the equivalence between the unitary condition and the re- 
quirement that the measurement of the ancilla extracts no information about the register 
system. {m\aij) = \{m\aij)\e'^^^^ ; the probability of the measurement does not distinguish 
between the register states and the term \ {m\aij) \ drops out and the effective unitary on the 
register pair of qubits is: 

U = diag(e^^°° , e^^°^ , e^^^° , e'"^'' ) (9) 

The result is also of the class A(a,0,0) and is thus also equivalent to a Control-unitary 
gate. Being diagonal in the computation basis, we can apply local unitary corrections to 
convert it into a gate of the form C — Rz{^). Given a two qubit unitary of the form 
diag(e*'^°°, 6*"^°^, 6*"^^°, e^'^ii) we can multiply it by local unitary gates 



= (iiag(e-^(«i+^i), e-^(^^+^2\ e-^^^2+^^\ 6-^(^2+^2^) 

(10) 

to give diag(e*^°°~^^^^^^\ g*0oi-(ai+&2)^ g«0io-(a2+^i)^ g*0ii-(tt2+&2) ^ 
We can choose ai, 6i, a2, 62 such that 

(/>oo-(ai+^i) = (11) 
(/>oi-(ai+^2) = (12) 
(/>io-(a2 + ^i) = (13) 
^ a2 + 62 = a2 + 61 - (ai + bi) + (ai + 62) = ^lo - ^00 + ^oi (14) 
the resulting gate ([9| must therefore be equivalent to the form 

U = diag(l, 1, 1, e^((^ii-^io)-(^oi-^oo) (15) 

We therefore will use ^ = d(j)i—S(j)o = (^11 — 0io) — (0oi — ^00) to characterise the entangling 
power of each gate. 



e 

ib 



e 
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4.2 A geometric picture of the unitary condition 

We have four final ancilla states \aij) corresponding to tiie basis states of the register. To 
allow a measurement basis {|m), \m_\_)} for which \ {m\aij) \ is constant Vz, j, these four points 
must all lie in the same plane and form a ring around a cap with |m) at the midpoint. For 
four \aij) given by the preparation of the ancilla state, the two ancilla-register couplings and 
the intermediate local unitary gate on the ancilla, |m) will be fixed and unique. The relative 
phases of (m|a^j) = \{m\aij)\e'^^^^ corresponds to the angles between the points on the ring 
(see figure [5|. 




Figure 5 For four states that have the same value \{m\aij)\, there are four points on the Bloch 
sphere that define a ring that encircle and thus define the state \m). \{m\aij)\'^ = cos^(^). 

4.2.1 A geometric picture of the entanglement condition 

After the first interaction but before the second the ancilla will be in one of two states \ai) 
corresponding to = z = 0, 1. The second interaction induces a unitary on the 
second register qubit that is given by {m\A\ai). For entanglement between the first and 
second register qubit, this unitary must be distinguishable by i. Therefore, each \ai) must 
have a distinct value of (cr^)^. In the geometric picture, two points given by that are on 
the same horizontal plane, before the second interaction, will produce four a^j on the same 
horizontal plane afterwards. However, the angle between states marked by different j in the 
pairs \aoj) and \aij) would be the same for each i. To be entangling Scjyi — (^0o must be 
non-zero. 



10) 10) 




ID ID 



Figure 6 All horizontal planes have a cap with a midpoint at the poles thus \m) = |0). An 
ancilla state in the computational basis corresponds to the same constant unitary enacted with each 
interaction and extracts no information from the first interaction to transmit in the intermediate 
stage. 
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From this geometric picture, we can show that the unitary and entangling conditions are 
fulfihed when the intermediate anciha states \ai) are restricted to the same vertical plane 
(see Appendix [A| and figure 12). 



4.3 Construction of the ancilla states 

With the intermediate states \ai) being of the same vertical plane but with some angle 
between them that we will label as 2/3, the final states \aij) will be 

(-^)'2^\0)+e^(-)-sinr^^i^) (16) 



COS 

We can construct the vertically split states \ai) and manipulate the angle 2/3 by the choice of 
preparation state and intermediate unitary Ua- Rotations around the z axis will form a gauge 
transformation since \m) is specified by {l^ij)} so we can, without a loss of generality, assume 
that the ancilla preparation state is on the x-z plane: \a) = cos(y)|0) + sin(Y)|l). After 
the first interaction this becomes \a) = cos(|-)|0) + e**^~^^'^"sin(|-)|l). Now the solid angle 
between the two resulting points is not given by 2a but by 2/3 where sin(2/3) = sin(6>')sin(2a). 
There will always exist a unitary that can rotate the two points such that they lie in the 
x-z plane: |a) = cos(^)|Q) +e^(-^)^'^sin(^)|l) ^ cos( ^~^~'^^'^ )|Q) +sin( ^~^~'^^'^ ). In the 
geometric picture, this is a rotation around the point where the great circle that connects 
the two points and the x-z plane cross, followed by any rotation around y of our choice so 
that ^ is a parameter under control (see figure [7|. The next interaction produces the points 
Eij which are, under this order of construction, dependent on the intermediate unitary and 
the ancilla preparation state. 



10) 10) 10) 10) 




ID ID ID ID 

Figure 7 A geometric representation of the method of constructing an entangling two qubit unit- 
ary. The first step corresponds to the ancilla preparation, the second and fourth to the interactions 
with register qubits and the third to unitary actions on the ancilla in between the interactions. 



4.4 An example of the relative entanglement powers of the Kraus 
operators 

As an example, take the intermediate state to have ^ = | so that the ja^) are vertically split 
about the |+) state. It is helpful to think of the intermediate state as Rx{-^)Rz{^2j3)\^) 
while the ancilla was prepared in |+). The effect of the preparation choice and the reduced 
solid angle are treated like an effective reduction in the interaction strength of the first 
interaction while the intermediate Ua = Rx{\)- The four states 



\a,j) = R,{{-iy2a)R,{-)R,{{-iy20)\+) 



(17) 



K. HalM Shah and D. K. L. Oi 



will be all symmetrically placed around |+) so we can also say \m) = |+) and measure in 
the {|+), |-)} basis. 



(m|a,,) = {MR,{{-iy2a)R, (|) i?,((-l)^2/3)|+) 



(18) 



R,(2a)R^(|)R,(2/3) 

^ (+|o^oo) 

(-l^^oo) 



1 f e-^^ -ie-'^ 



iB ^iA 



1 



(19) 



2^2 
1 

27i 



^ - ie-'^ - ie'^ + e'^) = ^(cos(A) - icos{B)) 

v2 



(20) 



le 



-iB I - iB 



ie*-" - e'^) = -^{-ism(A) - sin(B)) 



V2 



(21) 



where we define A = a-\-l3^B = a — (3. 
Each other element of the two qubit evolu- 
tion operator will just be a transformation of 
A and B (which thus fulfils the unitary con- 
dition) and because of the ± symmetries of 
sine and cosine we will be able to express the 
final Kraus operator and ^ in terms of just 
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^ Figure 8 By manipulating the effective coupling strength, /3, with the ancilla state preparation, 
the difference in $ for the two possible operator outputs can be adjusted. For a = f^, the difference 
between outputs is plotted with the solid (black) curve while the value of is plotted with the 
dashed (blue) curve (colour online). At /3 ~ 0.183, J$ = tt 



4.5 Repeat-until-success entangling gate generation (EGG) 

Without any optimisation of the ancilla initial state and measurement, the generation of 
entangling gates may, at the very least, be able to perform a random walk (a classical 
random walk, not to be confused with quantum walks) through the group of Control- 7^2(7)- 
With optimisation, we can apply a protocol to convert the output into a "success/fail" 
scenario with failure corresponding to generation of identity. 
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In this protocol EGG is enacted twice: once as described before and once again with 
local unitary changes that convert ^ — ^ — ^. The latter can be done either by local Pauli-X 
gate pre- and post-corrections on the register or, if there is access to the ancilla between 
interactions, by making a correction to the intermediate ancilla unitary to exchange the 
resulting intermediate states |ao) — \ai). The change in the second implementation enacts a 
change of the resulting gate C-Rzij) C-Rz{—^)- If these two output gates are (under local 
operator corrections) C-Rz{^q) and C-Rz{^i) with probabilities po and pi then after two 
two EGG the Control-Unitary is one of C-Rz{^q - <^i), C-Rz{^i - ^o) or I. If - = tt 
then we have performed probabilistic CZ generation with a success probability of 2poPi- 




C-R(Oi-Oa)=CZ 
=C-R(<I>o-Oi) 



H Figure 9 Since all enacted gates are of the same group, Control- i?2(<I>) gates, multiple gates can 
be easily combined. Random gate production is a Markovian process. By manipulating the relative 
values, the output can be limited to a finite probability tree, including a case equivalent to a single 
outcome "success/fail" gate suitable for a repeat-until-success method. 

For example, in figure |8) we see that we can lower x to match such conditions. In figure 
[To) the value of 2poPi is calculated, giving us a value oi p ^ 0.128. 
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■ Figure 10 The probability of achieving two opposing ancilla measurements, 2poPi, against /3, 
displayed against the resulting At the case where 5^ — tt, p ~ 0.128 (colour online). 

We provide this one adaptation as a single example of the general principle of using 
a multi-step protocol to adapt the behaviour of a Markov chain. Though an arbitrary 
interaction strength leads to the probabilistic generation of gates where the outputs are not 
local unitary correctable, they are of the same interaction class and thus can be adapted 
to the group of Control-Phase Rotations and generate group members. Specific physical 
parameter constraints and choice of target gates may lead to methods that involve different 
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schemes and different groups of gates such as C-R^i^)^ m^n G Z for a specific n or a 
continuous parameter group C-R^{'^). We provide a method that ahows one to enact a 
repeat-until success method similar to other probabihstic gate proposals for linear optical 
systems [2611271 128]. 

5 Initialisation-Measurement 

In ADQC, the ability to perform a projective measurement with an ancilla qubit was straight- 
forward and generalized measurements could be performed by introducing a second ancilla 
system. With a weak coupling, enacting a projective measurement is a less straightforward 
proposition requiring many steps and possible fidelity loss. We will provide a protocol for 
doing so in a naive adaptation of the Control- Z+Hadamard example, followed by discussion 
of potential future investigation and improvement. 

With a general register qubit = a|0) + and ancilla in initialised state 

the interaction E acting on \^)r\^)a would produce Q^|+)i?|0)A + /^|— A z basis 
measurement on the ancilla would replicate a z basis measurement on the register state 
before the interaction; this can also be used to initialise the register qubit with a result 
\j) producing a register qubit in the state Z^\^). Local unitary gate preparations on the 
register could be used to enable measurements or state initialisation in other abses. 

In our proposal, the interaction between register and ancilla can be weaker, enabling us 
to perform only non- projective measurements. Fortunately Oreskhov & Brun |29J show that 
weak measurements are universal and provides an explicit construction for how to use weak 
measurements to achieve projective and general measurements in a random walk, for two 
outcome measurements. This can be extended to higher dimensions [29 and it has been 
shown that one can build up many result general measurements from two output measure- 
ments [30 . Thus we expect that it should be possible to achieve at least an approximation 
of projective measurements using multiple successive non-project ive measurements from a 
single ancilla. We have a specific interest in being able to do so using our interaction of 
interest {Ha Hr).C-Tar- 

5.1 Initialisation and Measurement in logarithmic time with a fixed 
interaction 

At the heart of our proposal is the capability to perform a two qubit unitary that is equivalent 
up to local unitary gates to a Coniioi- Rz{0) gate. Consider |^) = a|0) and an ancilla 

in the |+) state: 



{HA^HR).C-R,{emR\^)A = H{a\{))^pco^[- |l))^|0)^ + /3sin - i7|l)^|0)A (24) 



Now if we measure the ancilla qubit in the {|0), |1)} eigenstate basis and make a local 
Hadamard correction on the register, the register qubit will be projected into the respect- 



ive states ^(q^|0) + /3cos(f )|1)) (with probability N = |ap -|- |/3pcos2(f)) and |1) (with 
probability |/3psin2(f )). 



Consider the effect of many such measurements in the case where every successive result 
on the ancilla is |0). The iterated effect over n steps is for the register qubit to be projected 
into the state ^=(a|0) + /3cos^(|)|l)) (where is the appropriate normalization factor). 
As n increases the state tends to |0) exponentially. The probability of achieving such a 
chain is Nn = |q^P + |/3pcos^^(|) which will tend to Thus we can replicate a projective 
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measurement in n steps up to a state fidelity error of [Scos^ ( f ) ^ cos^ ( | ) with the following 
procedure: 

Set a desired state fidelity error e. Calculate n G Z such that n > /n(e)//n(cos(|)). 
Prepare ancilla qubit in the |+) state, couple to register qubit with the interaction Ear^ 
measurement ancilla qubit in computational basis, apply H gate correction to register qubit. 
Repeat until measurement has been performed n times or until ancilla is measured as |1). 
Label a result of |0i0203...0n)A as |0)i^, any other result where the process was terminated 
by a \1)a as \1)r. 



|«p>- 



H 



■{}>>.■ 



H - H 



H 



T - H 



H 



H 



i3: 



T - H 



^ Figure 11 The circuit description of the iterative measurement protocol using the fundamental 
interaction Ear = {Ha Hr)C-T. The local unitary gate components of the interaction require 
a Hadamard gate correction; though not depicted, this is implemented using ancilla qubits as 
described in the previous section, leading to a total number of ancilla qubits and interactions of 2n. 

In this way, the measurement does not follow a random walk but occurs within a finite 
bound on the number of steps that is logarithmic with respect to the state fidelity error. 



^6 Avenues of further investigation 

Both the classes A(a,0,0) and A(aa;,a^,0) can be used to enact proportional-to-unitary 
Kraus operators however our work has mainly dealt with the former. Stochastic ADQC mo- 
tivates extending the method of enacting two qubit gates to the two parameter interaction 
class by providing means in which the /S.{ax^ocy^^) may be viable for single qubit gates. It 
also may provide a useful tool for examining the trade off between control or time. The 
random gate generation method we have described incoporates little control over the initial 
ancilla state or measurement basis. Potentially future work could consider employing feed- 
back to optimise the choice of ancilla state preparation and measurement basis in each step, 
resulting in a variable step size and guided random walk. The measurement we perform in 
each step is biased and we incorporate no feedback in our protocol save for the decision to 
continue or halt based on the previous measurement result. We have also only considered 
measurement on one qubit while the use of an ancilla provides us with a natural system 
for making higher dimensional measurements and generalised measurements. Combes & 
Jacobs [31] considered using feedback to speed up the purification of states via continuous 
measurements. A case where the measurement basis is kept unbiased to the state density 
operator through feedback is faster than measurement alone [32 . Much work has been done 
considering the case where we have a sequence of weak measurement with feedback enacted 
through unitary operations between steps in time for various goals and measures |3l] . 
Jacobs [35 raised the question of whether there is a trade off between speed of projection 
and loss of initial information. We propose the work of coming up with an unbiased basis 
feedback protocol for our specific interaction and non-continuous discrete protocol for future 
investigation. 
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l7 Conclusion and summary 

We have show how, by relaxing the requirement for stepwise determinism up to Pauli 
corrections, a broader class of interactions can be used to implement a form of ancilla 
driven quantum computation. This is achieved by gate approximation, repeat-until-success 
strategies, and generalised measurement. Single qubit unitary gates can still be implemen- 
ted deterministically while two qubit gates and measurement and state initialisation will 
be probabilistic and require the development of probabilistic protocols which we provide 
examples of. We provide a measurement process that is logarithmic to the bound on state 
fidelity error and reduces to a deterministic approximation of projective measurements, with 
the potential for further speed up to be examined in future research. The conditions for 
the generation of a two qubit entangling unitary gate are described in geometric terms. We 
demonstrate how choice of ancilla state preparation and local unitaries enable us to manip- 
ulate the entangling strength of the output gates. This allows us to convert the result into 
a "repeat-until-success" gate. Further research may reveal other protocols based around 
adaptations of the entangling gate power. Additionally, probabilistic protocols for enacting 
specific unitary gates may also aid in expanding possible interactions to include the class 
A(a^,a^,0). 

We would like to acknowledge discussions with Dan Browne, Erika Andersson and Shojun 
Nakayama, and comments on the draft manuscript by Janet Anders and Elham Kashefi. 
DKLO and KHS are supported by the Quantum Information Scotland (QUISCO) network. 
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Appendices 



A Restrictions of the parameters of the ancilla found in a geometric 
proof 

Each point can be represented by the spherical coordinates that describe the state cik = 
{^k^(t^k)i k = 1, 2, 3, 4. For out notation the values of k will correspond to the values expressed 
by ij in binary. 

Each point can be expressed in Cartesian coordinates by the relationships: 

\ak.x\ = sin(6>fe)cos(0fe) (25) 
\dk-y\ = sin(6>fe)sin((/)/e) (26) 
\ak.z\ =cos(6>fc) (27) 

Three points alone can always be found to be on the same plane. We will define a plane 
from three points and then find the expression for the distance from the fourth point to that 
plane. Thus we will find the conditions for the fourth point to be in the same plane as the 
other three. The equation for a plane defined by three points is 

a.x^b.y^c.z^d = (28) 
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(29) 



There is a freedom of choice of d so it can he sei d = D. The distance of point 4 to the 
plane is given by 



distance ■ 



|a.X4 + 6.7/4 + + d\ 
+ 62 + c2 



(30) 



As we are only interested in the case where the distance is zero, we can ignore the normal- 
isation factor and simply examine 



distance' = |a.X4 + h.y/^ + C.Z4 + d\ 



(31) 



We now make use of some restrictions on the formation of these four points. Before the 
second interaction, there are intermediate ancilla states corresponding to the computational 
basis of the first register qubit. The final four points are constructed from these two points 
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by the second interaction which imparts a rotation around the z axis by ±2a with ± cor- 
responding to the second qubit being in the computational basis state j = 0,1. This 
means that we can set 



(32) 



and 



4a 



(33) 



Using ( 32 ) , the distance between the fourth point and the plane can be simplified to 



2(cos(i92)-cos(i94)) 



cos 



Z>3 + 04 
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sin(^2)sin(^4)sin(- 



(34) 



This generates several possible conditions for the fourth point to lie in the plane, some 
more trivial than others. If cos(6>2) = cos(6>4) then all four points must lie on the same 
horizontal plane which means that there is no entangling power. sin(6>2) = and sin(^4) = 
would mean that there are only three distinct points with one being at the pole i.e. the |0) 
state. Due to the construction of these four points it is not possible for ^3 — ^4 = to be 
true. The final condition is that 



cos 



h + 04 
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If (33) is then substituted in, this becomes 
203 +4a^ 
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(35) 

(36) 
(37) 



Since a is non-zero, this requires 0i 
vertical plane. 



03 - 2a) 

F nTT for n G Z i.e. the two points are in the same 




" Figure 12 A 2d projection of the construction of four point on the Bloch sphere. By restricting 
the points to only one of two polar and one of two azimuthal angles, the vectors that connect two 
points of the same polar angle will be parallel. This guarantees that all four points lie on the same 
plane. 



B Graphs of stochastic ADQC gate count and the exponential 
distribution 
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Figure 13 Probability distribution of required gate count to achieve target Ut = Rx{^)- a) Use 
of a single parameter interaction in a 20 bin histogram, b) Use of a two parameter interaction in a 
20 bin histogram, c) Use of a single parameter interaction in a 100 bin histogram, d) Use of a two 
parameter interaction in a 100 bin histogram; note the large aberration in the first division. The 
probability distribution corresponding to the exponential distribution parametrised by the mean of 
the results is displayed by the solid red curve (colour online). 
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P Figure 14 The natural logarithm of the bin counts of the number of gates required to achieve 
target Ut = i?x(f )• By taking the natural logarithm the exponential distribution model forms a 
linear curve which fits the points generated by the simulation. Noise distrupts the linear behaviour 
for very low counts of high gate number. Figure a matches the one parameter interaction, figure b 
matches the two parameter interaction. 



